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Plug-in electric vehicles (PEVs) are expected to balance the fluctuation of renewable energy sources 
(RES). To investigate the contribution of PEVs, the availability of mobile battery storage and the control 
mechanism for load management are crucial. This study therefore combined the following: a stochastic 
model to determine mobility behavior, an optimization model to minimize vehicle charging costs and 
an agent-based electricity market equilibrium model to estimate variable electricity prices. The variable 
electricity prices are calculated based on marginal generation costs. Hence, because of the merit order 
effect, the electricity prices provide incentives to consume electricity when the supply of renewable 
generation is high. Depending on the price signals and mobility behavior, PEVs calculate a cost minimizing 
charging schedule and therefore balance the fluctuation of RES. The analysis shows that it is possible to 
limit the peak load using the applied control mechanism. The contribution of PEVs to improving the 
integration of intermittent renewable power generation into the grid depends on the characteristic of 
the RES generation profile. For the German 2030 scenario used here, the negative residual load was 
reduced by 15-22% and the additional consumption of negative residual load was between 34 and 52%. 

© 2012 Elsevier Ltd. All rights reserved. 
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1. Introduction distributed load and generation units, we distinguish between 


direct control and indirect control. Direct control or centralized 


Integrating high shares of fluctuating power generation into the 
electricity system requires flexible power plants, storage, power 
distribution via a reliable power grid and greater flexibility on the 
demand side. In terms of demand response [1] or the control of 
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optimal charging implies that a service provider can shut down or 
reduce loads and control decentralized generation units directly. 
Examples include the direct load control of residential water 
heaters [2] and air conditioning loads in California, or of virtual 
power plants such as the “Dezentrales Energie Management 
System” (DEMS) of the German Siemens AG. The advantages of 
such direct control are prompt and predictable reactions to control 
signals. Drawbacks arise from reduced consumer acceptance in 
the case of controlling loads in private homes or vehicles and the 
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communication and optimization efforts involved in controlling a 
large number of small storage or generation devices with varying 
consumer needs and providing mobility as primary purpose of use. 
Indirect control uses price signals to manage loads or generation 
units. In this case, the service provider sends price signals and 
the consumer (or an automatically controlled device programmed 
by the consumer) decides to either reduce or shift the load when 
the price is high, or pay the higher price. In this case, consumer 
acceptance should be higher than is the case for direct control 
[3]. Disadvantages arise from the possibility of avalanche effects 
or simultaneous reactions [4,5] and the necessity to predict the 
reaction of consumers to different price signals, which yields 
the possibility of forecasting errors. However, since consumer 
acceptance seems to be crucial for the feasibility of managing 
mobility-related systems, an indirect energy management system 
is considered the most promising option to control plug-in electric 
vehicles (PEVs). Hourly prices in combination with smart devices 
(devices such as thermostats or electric vehicles that optimize 
their demand depending on a control signal) are an adequate 
tool to involve consumers in the electricity markets [6,7]. The 
contribution of PEVs to improving the integration of intermittent 
renewable energy sources (RES) into the grid depends on technical 
issues such as storage capacity, grid connection power, and driving 
behavior, which together define the energy available for load 
shifting, as well as social and economic aspects which influence 
the incentive for consumers to participate in the load-shifting 
program. In this paper, we focus on modeling technical issues and 
driving behavior (see Section 2.3). Several recent studies discuss 
the impact of PEVs on power systems but do not consider the 
impact of RES on load shifting strategies [8-10], the dynamic 
mobility behavior and/or load shifting mechanisms [11-12]. 

For the analysis presented here, the agent-based electricity 
market equilibrium model PowerACE (see Section 2) is used to 
determine dynamic prices for a German 2030 scenario with a high 
share of intermittent generation (for the assumptions, see Section 
3). The intermittency ofthe renewable generation affects the price 
signals (see Section 4) used for the distributed optimization of PEVs 
charging schedules. Hence, PEVs preferentially consume electric- 
ity when the supply of intermittent RES is high. Avalanche effects 
that can occur because of the similar optimization results of many 
devices are discussed in Section 5. In a case study for Germany, 
the effect of demand-side management with PEVs on the residual 
electricity load is evaluated and their contribution to integrating 
intermittent RES in the grid is shown (Section 6). 

The present work is closely related to the Volkswagen and E.ON! 
field test “Flottenversuch Elektromobilitat” funded by the German 
Federal Ministry for the Environment, Nature Conservation and 
Nuclear Safety (BMU). In the field test, price signals generated with 
the PowerACE modelare sent to Volkswagen plug-in hybrid electric 
vehicles (PHEVs). Further, the same algorithms applied to optimize 
the charging behavior of the vehicles in dependence on price signals 
are also used in the simulation. In this way, consumer and real-life 
driving aspects will be able to be addressed in the future. 


2. Simulation approach 


The effects of PEVs on the electricity system are investigated 
using the market equilibrium model PowerACE in combination 
with indirect energy management and stochastic modeling of driv- 
ing behavior. On the demand side, static and dynamic electricity 
demand as well as grid losses are considered on an hourly basis 
over one year of simulation. The dynamic demand is related to 


1 E.ON is utility enterprise. 


distributed devices that are able to optimize their market behav- 
ior, suchas night storage heaters, electric vehicles or other shiftable 
loads. 

On the supply side, RES, conventional power plants and storage 
technologies are considered. Details of the PowerACE model can be 
retrieved from [13] and an overview of the model is given in Fig. 1. 


2.1. Demand-side management agent 


Electric vehicles are implemented as part of a dynamic demand 
agent. The dynamic demand agent is structured in pools, groups 
and devices. A pool bundles several or at least one group and the 
group bundles several or at least one device.? In this study, a device 
represents a PEV. The group level is used for regionalization and 
grid restrictions. The device level represents all vehicle-specific 
information, including driving behavior and price-dependent opti- 
mization (see Section 2.2). Vehicle pools administer data for the 
groups and devices belonging to the specific pool. Each pool acts 
as an independent agent and can place demand and supply bids in 
the spot market. A pool adds up the expected operation og, of all 
the relevant devices in the different groups G. To keep the model as 
simple as possible, the operation o of a device din group g is known 
at the pool level when the bid is placed in the spot market.? The 
operation Op of pool p is calculated using Eq. (1): 


Dı Dz Dg G Dg 
Oplt)= X oag, (t)+ X` oag lt) + XO oaclt) = XX 045,¢-(1) 
dz1=1 dya= dgn=1 g=1 dg=1 


As a guidance signal for the devices participating in a pool, each pool 
performs a price forecast. The price forecast and the supply bids 
which result in the market clearing price are calculated according 
to the marginal costs. The marginal electricity costs Cmarginal Consist 
of the fuel price Prueı to produce 1 MWh of electricity and the costs 
for the resulting CO, emissions. The CO, costs are given by the CO 
coefficient, qco,, which defines the CO, emitted when transforming 
a primary energy carrier, and the CO; price pco, . Dividing the prices 
for fuel and CO; by the efficiency n of a specific power plant gives 
the marginal costs: 


1 1 
Cmarginal = 7 * Pfuel + 7 * dco, * PCo,- (2) 
To calculate the supply bid (bidsuppiy price/capacity) of a power 
plant with a specific capacity, a markup factor mup is added to the 
marginal costs: 


. Cmarginal + Mup 
bidsupply = ~ Capacity (3) 


The markup factor represents the margin or a markup to cover total 
costs. In the work presented here, base load power plants (on the 
left-hand side of the merit order in Fig. 2) in some cases place bids 
below the marginal costs to avoid start-up operations (negative 
markup), and peak power plants (on the right-hand side of the merit 
order in Fig. 2) place bids with positive markup to cover the full 
costs of the power plant. Because of the very low costs of operation 
for intermittent RES, marginal electricity costs of zero are assumed. 
Hence, RES supply with priority and reduce the power required 
from controllable power plants (Section 4). 

The 2030 market clearing prices (PowerACE) for one year 
(8760h) and the merit order of all power plants available in the 
2030 scenario (see Section 3) are shown in Fig. 2. The figure also 


2 The structure was selected to model a general demand-side management agent 
which could not only include PEVs but also devices with thermal storage such as heat 
pumps and freezers or DSM devices without storage. Here, only PEVs are included. 

3 The pool places a price-independent bid. Hence, the bid must be always cleared. 
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Fig. 2. Merit order and clearing price bids in a German 2030 scenario. Simulated with the PowerACE model; fuel and CO3 prices according to [19] Scenario A “deutlich”; 
installed convention generation capacities based on own estimations. Installed capacity: oil 0.7 GW; gas turbines 16.4GW; combined gas and steam 10.2 (n=45-59%) and 


25.6 GW (n=60-65%); coal 8.7 GW, lignite 9.2 GW; waste 0.9 GW. 


gives a price forecast function of the marginal electricity costs over 
the residual load (RS) which is used to determine the forecast price 
ofa pool. To calculate the residual load RS (Eq. (4)), perfect foresight 
for the inelastic electricity demand Dynelastic, the exchange balance 
Ex and the generation of renewable energy sources RES are used: 


(4) 


The expected operation of all pools Op is calculated individually for 
each pool one after the other according to Eq. (5): 


p-1 
P 
Op p(t) = SM). 
n=1 


The operation of the pools in which devices have already opti- 
mized their demand (pool p) is known and scaled up to estimate the 
operation of all pools. Using the residual load, the pool performs a 
two-day forecast for the expected prices Pfforecast USING Eq. (6): 


RSp(£) = Dinetastic(t) + EX(£) + Op(t) — RES(t). 


(5) 


Pl forecast,p(RS) = 0.0071 - RS? — 0.4723 - RS? + 10.813 RS + 10 
(6) 


for RS values greater than zero. For a RS equal zero or negative a 
linear correlation‘ is used. 

The assumption that the generation of intermittent RES and the 
elastic demand is known is obviously not realistic for an actual elec- 
tricity market. This method is used because the focus of this analysis 
is the indirect control of PEVs and their ability to integrate inter- 
mittent RES. A study discussing possible market organizations is 
presented in [14]. 


2.2. Vehicle-based optimization of the load profile 


The used simulation structure? considers the superior market 
situation represented by the pool level and the local distribu- 
tion network situation represented by the group level. Devices are 
affected by both levels due to different tariff components. 


4 pr=RS+10. 
5 Simulation structure, a pool contains groups and a group contains devices or 
PEVs, respectively. 
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At the group level, a variable grid fee is added to the price com- 
ponent of the pool level to avoid overloading in the distributed 
network (see Section 5). Therefore, a simple quadratic relation 
between the variable grid fee prgridfee,g and transformer utilization 


Utransformer,g 15 used: 
2 
Pl gridfee,g(t) = A * UT, ancformer,g + ©: (7) 


The constant parameters a and c are selected, such that the total 
grid fee is constant for one day at a transformer. Hence, the sum of 
the constant grid fee multiplied by the power at the transformer for 
all time steps t of a respective day is the same as using the variable 
grid fee assuming that the vehicles charge after their last journey 
of the day. 

The transformer utilization is calculated by dividing the power 
at the transformer in time step t by the nominal transformer power 
Wrransformer- The power at the transformer is calculated by adding 
the standard load profile W(t) for Germany [15] and the operation 
Oq g(t) of the PEVs already processed: 


Wit) + 04 g(t) 


Wrransformer 


UTransformer,g(E) = . (8) 
After each optimization of a device, the operation 0q g(t) of the PEVs 
is known and the expected utilization of the transformer changes. 
As for the pool price forecast, this is a theoretical approach to avoid 
avalanche effects (see Section 5). Today, equal treatment of electric- 
ity consumers is required and therefore customer-specific tariffs 
are not allowed in Germany. 

The resolution of the grid price and the operation of PEVs at grid 
level are 15 min blocks of time or 96 time steps per day. The devices 
optimize their operation depending on a price signal to minimize 
costs. The price for a specific device pry comprises the pool forecast 
price Pfforecast,p (See Eq. (9)), the variable grid fee Pgriafee,g and the 
sales tax Itax: 


Pra(t) = (Pfrorecast,p(E) + Pfridfee,g(E)) x (1 + tax). (9) 


The sales tax used do not influence the results in the simulation. 
This value is implemented because the calculated tariff is also used 
in a field test [16] and the optimization is done from a consumer 
point of view. In this case, the sales tax widens the price spread and 
therefore increases the consumer incentives. 

The local optimization uses a graph search algorithm to mini- 
mize the charging costs [17]. The algorithm calculates a schedule 
for charging and discharging the vehicle battery in the optimization 
time period. The algorithm is applied when the vehicles are plugged 
into the grid after each trip. In addition to the price prg, the usable 
battery capacity, the state of charge (SOC), the beginning and end of 
the optimization, the grid connection power and battery degrada- 
tion all enter the schedule calculation. After a new optimization, the 
previous schedule is invalid. A detailed description of the applied 
algorithm is given in [18]. 

The optimization time period depends on the SOC when return- 
ing to the grid and follows the values in Table 1. 


2.3. Driving behavior 


A stochastic approach is used to estimate driving behavior. 
The data is taken from the travel survey “Mobility in Germany” 
[19]. In order to find user segments that are suitable for electric 
driving, the raw data is filtered [20]. Early adopters of PEVs can be 
characterized by a high annual driving mileage due to economic 
reasons. Compared to conventional vehicles the investments in 
PEVs are higher but the operating costs are lower. We consider 
technical aspects as well as economic aspects. Possible PHEV and 
BEV users are required to have an assigned parking space. Addi- 
tional requirements for BEVs include the availability of more than 


one vehicle per household and regular job-related trips shorter 
than 90 km. 

Modeling driving behavior can be simplified to three probability 
parameters which are defined as follows. 

The probability to travel with the vehicle on a certain day 
Protrayeı (day): 


Maga, 
Prorravei (day) = x x "Travel;(day). (10) 


Mday 
where Travelaay is a Boolean value (true, false) indicating whether 
the respondent is driving on day and mgay represents the sample 
size on a given day of the week. We found distinct patterns of trav- 
eling behavior on Sun, Sat, Mon, Fri and other weekdays. Weekdays 
are defined as Tuesday, Wednesday and Thursday and merged into 
one data set because the driving behavior on these days was found 
to be very similar. 

The probability for starting a trip Prostar: on a specific day and 
time slot is given by: 


r I" nday, t). (11) 


Naay represents all samples with a trip on a specific day and n(day, 
t) a trip started at a specific time t. t is a parameter out of 0-95, 
or a 15 min time resolution during a day, respectively. The average 
trips! per day, while Pro;,aveı = true are multiplied with Probstart to 
account for multiple trips per day.® 

Range and duration are correlated. The different range values 
are Classified in k:= {0, . . ., 20}. The probability Prorange at a specific 
day for a range class k is: 


PrOstart(day, t) = Na 
ay 


x) "ntday, k). (12) 


Prorange(day, range € k) = i 
i= 


Nday 
The class specification is given in Table 8 in Appendix A. 
For the data used, a linear correlation was found between range 
r and duration du, which is shown in Fig. 3. 
The duration du(r) is given in Eq. (13): 


du(r) = 0.6837r +5. R? = 0.92 (13) 


Including aspects related to the infrastructure in the simulation 
model requires an additional probability value: 


Nday 
x > n(day, t, k). (14) 


Pro,..(day, t, location € k) = 


Nday 


PrOjocation represents the probability trip ending at k out of 0... 4 
locations. In the work presented here, infrastructure is not consid- 
ered. We assume that PEVs are plugged-in after each travel trip and 
that the necessary infrastructure is available. 


3. Assumptions 


In order to analyze the effect of fluctuating renewable energy 
generation from wind power and photovoltaic, as well as the 
contribution of PEVs towards balancing this RES, we constructed 
a scenario for 2030. The hourly characteristic of the demand 
in Germany is taken from [21]. For 2030, a yearly demand of 
502.1 TWh is assumed [22]. To indicate the supply side, the merit 
order of power plants is generated using primary energy and CO, 
prices from [23, Scenario A “deutlich”]. Based on the current Ger- 
man power plant mix, all power plants reaching the end of their life 
cycle life by 2030 are assumed to be replaced by gas turbine power 
stations and combined-cycle plants. The capacity of intermittent 
renewable energy sources (RES) for the 2030 scenario is taken from 


6 Average trips per day Mon=3.950, Weekday =3.965, Fri=4.211, Sat=5.597, 
Sun = 2.810. 
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Table 1 
Optimization time to recharge the battery storage. 
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Optimization State of charge (SOC) SOC <0.3 SOC <0.7 SOC = 0.7 Unit (1 = 15 min) 
After a trip Storage <5 kWh 8 32 60 Time steps 
Storage > 5kWh 16 56 60 Time steps 
Every day at 22 o’ clock - 38 38 38 Time steps 
Table 2 
Used structure of pools, groups and devices?. 
Pool Group Dev1 PHEV27 Dev2 PHEH57 Dev3 City-BEV Dev4 BEV Sum Grid power Total storage Transformer capacity 
[1000 units] [GW] [GWh] [kVA] 
1 1 129 220 41 10 400 1.80 4.14 AkVA/Dev + Household 
1 2 130 219 41 10 400 1.80 4.13 AkVA/Dev + Household 
1-9 Sum 2331 3951 738 180 7200 32.47 74.37 - 
10 1 259 439 82 20 800 3.61 8.26 AkVA/2x Dev + Houshold 
10-15 Sum 1554 2634 492 132 4812 21.74 49.94 - 
1-15 Sum 3885 6585 1230 300 12,000 54.12 123.95 - 
a Source: Dominanz-Szenario [20]. 
600 
+ Range and average duration (MiD) . . 

= 500 p 

£ — Linear correlation between range and duration _- 6 

E w _ d z 

a = + 

3 a wos . 
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Fig. 3. Correlation between the average duration of a trip and the range of a trip. The used data is filtered according to [20] and uses raw data from [19]. 


[23,Scenario A]. The installed capacity of onshore wind, offshore 
wind and photovoltaic is 37.8 GW, 25 GW and 63 GW, respectively. 
The hourly characteristic of intermittent renewable energy sources 
for electricity production is taken from [21] for wind onshore and 
[24] for wind offshore and photovoltaic. The imports and exports 
of electricity and storage technologies like hydro pump storage are 
not taken into account. 

The demand of BEVs and PHEVs is modeled (see Section 2) using 
the probabilities describing the driving behavior of the mobility 
survey MiD 2002, which was aggregated according to [20]. The 
probabilities used are available in Appendix A. 

The penetration of PEVs is taken from [25]. In the 2030 scenario, 
12 million PEVs will be on the roads in Germany. Table 2 shows the 
allocation of the PEVs in vehicle pools. The pools 1-9 are composed 
of two groups and pools 10-15 are composed of one group which 
results in a different transformer capacity available per vehicle. In 


Table 3 
Assumptions for the different types of plug-in electric vehicles (PEVs). 


Dev Type Storage usable Grid connection Equivalent energy 
energy [kWh] power [kW] use [kWh/km] 

1 PHEV(25) 4.5 4 0.18 

2  PHEV(57) 12 4 0.21 

3 _ City-BEV(100) 15 8 0.15 

4 _ BEV(167) 30 8 0.18 


PHEV: plug-in hybrid electric vehicle and BEV: battery electric vehicle. Value in 
brakes is the theoretical driving range in km 


this model, the operation of the vehicles is scaled-up by a factor 
of 1000. In total, 12 thousand PEVs are modeled, representing 12 
million PEVs. 

The different PEV types used are shown in Table 3. The values 
represent an average of the specific vehicles. 

Inthe quadratic equation to estimate the grid-related price com- 
ponent (see Section 2.2), a is assumed to be 30 and c is 0.5. 


4. Characterization of the merit order effect 


An important aim of the work currently being conducted is 
to investigate how PEVs can be used to improve the integration 
of intermittent RES and balance the fluctuation of these power 
sources, mainly wind power and photovoltaic. The effect of inter- 
mittent RES on electricity prices, known as the merit order effect, 
plays an important part in the applied approach [26]. The indirect 
control of PEVs uses a tariff that depends on marginal electricity 
costs and therefore on the residual load. 

Fig. 4 shows the merit order of the power plants available in 
Germany in 2008 (left curve in Fig. 4). Taking the same marginal 
costs and power plants, but including the installed capacity of wind 
power and photovoltaic shifts the merit order to the right (right 
curve in Fig. 4). Itis unlikely that the total installed capacity of inter- 
mittent RES is available simultaneously, but it can be concluded that 
these intermittent RES do reduce the marginal costs and therefore 
the clearing prices in the electricity market at least to some extent 
[27]. 
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Fig. 4. Merit order of the power plants in Germany 2008 with and without the installed capacity of wind and photovoltaic. 
Source: own calculation using data from [23] for the installed capacity of wind and photovoltaic. 
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Fig. 5. Merit order of the power plants in a 2030 scenario for Germany with and without 50% of the installed capacity of wind and photovoltaic. 
Source: own calculation using data from [23] for the installed capacity of wind and photovoltaic. 


For the 2030 scenario (see Section 3), the effect on the merit 
order due to photovoltaic and wind generation increases (see 
Fig. 5). We assume a situation in which energy generated from 
photovoltaic and wind makes up 50% of the installed capacity. This 
situation will result in marginal costs of zero for 6(0 GWh demand.’ 
If no intermittent RES supply is available in this demand situation, 
the marginal costs would be 120 € per MWh. 

This example shows how intermittent RES affect the marginal 
costs of electricity and that future electricity markets will have 
more volatile prices. It is very likely that higher shares of intermit- 
tent RES will also increase the base peak spread of electricity prices. 
This is because of an expected increase in fuel and CO, prices and 
a tendency for bids to include a higher markup factor to cover the 
full costs of conventional power plants. 


5. Avalanche effects 


Indirect control of distributed generation and flexible loads with 
variable tariffs has one major drawback. If an automated optimiza- 
tion of devices is used, a price signal which is valid for all devices 
can cause demand or generation peaks, known as avalanche effects 
[4,5,28]. These effects are induced by an optimum to consume or 
feed-back electricity in one specific time period which has the low- 
est or highest prices for electricity. Assuming that all devices have 
the same degree of freedom (actual and favored state of charge at 


7 The German simultaneous hourly annual peak load in 2008 was 76.8 GW. 


the end of the optimization time, grid connection time, and soon) in 
a specific time period for load-shifting or generation will reinforce 
this effect. In the case of electric vehicles, however, the degree of 
freedom is affected by driving behavior and the willingness of the 
users to shift loads and to plug-in the vehicle.® 

Fig. 6 shows the optimization result of PEVs for three PEV pools 
using the same price forecast (no grid fee is assumed). Driving 
behavior affects the amount of energy required to recharge the 
battery and the time steps in which the PEVs are connected to the 
grid. The demand of all PEVs collected in the pools 1-3 is very sim- 
ilar. Only pool 1 shows slightly different results. The diffusion of 
the demand results from the limited grid connection power? and 
differences in driving behavior. 

Using individual price forecasts for the different pools results in 
a greater diffusion of the PEV pool demand (see Fig. 7).! 

The avalanche effects are reduced and PEV demand is better dis- 
tributed to fill load valleys (see Fig. 9 in Section 6). In real electricity 
markets better allocation is achieved using price elastic demand 
bids.!! 


8 The willingness of the consumer to react to incentives is not part of this paper. 
We excluded this aspect in order to examine the effects of PEVs without the neces- 
sary assumptions about consumer incentives. 

9 If it is not possible to completely recharge the battery in the time step with the 
lowest price, the next time step is used. 

10 Perfect forecast is assumed. Hence the reaction of other pools is known (feed- 
back) and/or estimated (see Section 2). 
11 [10] shows that a Nash equilibrium in case of PEVs is possible. 
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Fig. 6. Demand of three PEV pools with the same price signal. 


0.3 
——Price forecast (Pool 1) 


0.25 


Load [MW] 
9 
i 


0.05 


— Load Pool 1 


— Load Pool 2 


— Load Pool 3 


Price forecast (PowerACE) [€/MWh] 


1 5 9 131721 1 5 9 131721 1 5 9 131721 1 5 9 131721 1 5 9 131721 1 5 9131721159 
Time [h/day] 


Fig. 7. Demand of three PEV pools with different price forecast. 


Avoiding peaks in the total system load can still result in load 
peaks in distribution electricity networks. A strong increase in 
wind generation, for example, can necessitate the concentration 
of additional demand within a short time period. At the group or 
distribution network level, therefore, a variable grid fee is used to 
guarantee that the transformer is not overloaded. This variable grid 
fee is specific to each device (see Section 2). Due to the assumed per- 
fect demand forecast of other devices and electricity users within 
the distribution network, the grid fee increases after each demand 
optimization of a device (feedback loop). Hence, the optimal load- 
ing time period shifts and the utilization of the transformer is 


F(x)=ax?+bx+c 
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evened out. Depending on the situation at the system level and 
in the local network, price components at the pool and the group 
level can provide incentives for the same time periods or counteract 
each other. 

Fig. 8 shows transformer utilization with and without PEVs for 
one group (we assume that a 400 kVA transformer hosts 100 house- 
holds and 100 PEVs). The function of the dynamic grid fee over the 
transformer utilization is also shown. 

Transformer utilization is not critical with the variable grid fee 
used for the 2030 scenario. An increase from about 25% to 35% on 
15 min basis can be observed. 
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Fig. 8. Transformer utilization with and without plug-in electric vehicles (PEVs) and the function between the dynamic grid fee and transformer utilization. 


D. Dallinger, M. Wietschel / Renewable and Sustainable Energy Reviews 16 (2012) 3370-3382 


100,000 
@ PEVs last trip 


60,000 


40,000 


Load [MW] 


20,000 


3377 


Total system load ™PEVsDSM Residual load 


Fig. 9. Load of plug-in electric vehicles (PEVs) charging after the last trip and using demand-side management (DSM). First three weeks of the year 2030 scenario. Renewable 


energy fluctuation based on the weather year 2008 [24]. 
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Fig. 10. Sorted total system load and sorted residual load 2010 as well as sorted residual load 2030. See Section 3 for the assumptions. 


6. Results 


The model results based on the assumptions outlined in Section 
3 are presented in the following section. Unlike other publications 
(e.g. [8,9]), the PEVs’ load follows the residual system load and does 
not simply shift demand to the night-time hours (see Fig. 9). The 
intermittency of renewable energies makes it necessary to take this 
advanced view of base and peak load. Further, the load shifting fol- 
lows a dynamic pricing approach with a distributed vehicle-based 
optimization which considers mobility behavior, local grid restric- 
tions and consumer incentives. 

Comparing the demand-side management (DSM) with simple 
charging after the last trip in Fig. 9 shows a reduced peak load. 
PEVs preferentially consume energy when the residual load is low 
or even negative. 

In the assumed scenario, the installed capacity of wind and pho- 
tovoltaic rises from 46 GW in 2010 to 125 GW in 2030 (see [23], 
pp. 186-187). Comparing the residual load in 2010 with the resid- 
ual load in the 2030 scenario (see Fig. 10) illustrates the switch to 
electricity systems with higher shares of intermittent generation 
from renewable energy sources.!? For 2030, we expect 635 h with 
a negative residual load (in total 5.36 TWh), which means that the 
generation of intermittent photovoltaic and wind power plants is 
higher than the assumed demand. The minimum of the residual 
load drops from 10.4 GW in 2010 to —43.7 GW in 2030. 

In order to analyze the PEVs’ contribution towards improving 
the integration of intermittent RES, we focus on three key values: 


12 In the simulation presented here we did not consider pumped storage or trans- 
border flows. 


the change in the minimum residual load, the percentage of neg- 
ative residual load that can be consumed using load shifting, and 
the residual load change (RSt) — RS(r- 1)). The characteristic of the 
RES intermittency and the energy produced influences the PEVs’ 
ability to balance RES. To account for this we analyze different 
weather years with specific wind on- and offshore as well as pho- 
tovoltaic data for 2006-2008 [21,24]. To compare the different RES 
characteristics, we additionally scale the 2006 and 2007 data to 
the RES energy production (325 TWh) in 2008 (see Tables 4 and 5 
“scaling”). 

Fig. 11 shows the effect of PEVs on the residual load for the two 
charging scenarios with weather characteristics for 2008. Com- 
pared to uncontrolled charging after the last trip, demand-side 
management (DSM) can prevent an increase in the peak load (for 
all the weather years analyzed). 

The minimum of the negative residual load is reduced by 8.0 GW 
to —35.7 GW in the case of DSM (see Fig. 12). 2.82 TWh or 52.6% of 
the negative residual load can be consumed using load shifting. 

Considering the different weather years in Table 4 shows a peak 
reduction in the negative residual load of between 15 and 22% or 
about 6-8 GW. The negative residual load usable with DSM (=(neg. 
RS — neg. RS DSM)/neg. RS) varies between 34 and 52%. This corre- 
sponds to 14.3% (2.56 TWh) and 20.4% (3.66 TWh), respectively, of 
the 17.9TWh consumed in total by the PEV fleet. 

This shows that the characteristic of intermittent RES strongly 
influences minimum peak reduction and negative RS able to be used 
by DSM. A more balanced RES generation with a higher availability 
on each day of the year (e.g. offshore wind) or short recurring pho- 
tovoltaic generation peaks are better able to be used than extreme 
weather conditions with particularly strong or no wind periods. 
That is because, in the case of DSM, only the electricity used for 
driving (on average 49 GWh per day) can be used for load shifting. 
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Table 4 

Results of the negative residual load (RS) peak reduction and usable negative residual load with different renewable energy intermittency. 

Unit 2006 2006 scaling 2007 2007 scaling 2008 

RES generation Twh 306.97 325.20 336.87 325.20 325.10 
Min. RS power GW -37.19 —43.99 —40.75 -37.35 —43.74 
Min. RS power DSM GW -30.72 37.05 32.01 28.90 35.74 
Min. RS reduction DSM 17.42% 15.78% 21.45% 22.61% 18.29% 
Negativ RS Twh 6.16 9.87 9.25 6.83 5.36 
Negativ RS DSM Twh 3.60 6.44 5.59 3.89 2.54 
Negative RS usable DSM 41.55% 34.77% 39.61% 43.09% 52.58% 


Table 5 
Results of the residual load (RS) change versus RS change including plug-in electric vehicles (PEVs) using demand-side management (DSM) with different renewable energy 
intermittency. 
Unit 2006 2006 scaling 2007 2007 scaling 2008 
Residual load (RS) 
Load change Min. GW —15.39 —16.16 —14.64 —14.15 —14.41 
Load change Max. GW 18.88 19.77 16.33 15.82 22.13 
RS + PEVs (DSM) 
Load change Min. GW 15.71 16.76 12.88 11.98 — 14.45 
Load change Max. GW 18.73 19.60 16.05 15.77 17.90 
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Fig. 11. Effect of plug-in hybrid electric vehicles (PHEVs) on the residual load when charging after the last trip and using demand-side management (DSM). See Section 3 for 


the assumptions. 


In general, mobility behavior restricts load shifting capability, so 
PEVs are suitable for short-time storage (1-2 days). 

The hourly load change of the residual load with and without 
price-controlled PEVs is shown in Fig. 13. The effect of different 
weather years on the residual load change is relatively small 
(represented by the width of the lines in Fig. 13). Further, for the 
RES generation and load curve considered, the average positive 
delta (RS) > RS¢t_ 1) on the right side of Fig. 13) is higher than the 
average negative delta (RS, < RS;;_ 1) on the left side of Fig. 13). As 


a result, a rapid generation increase is required more often to bal- 
ance the system. Price-controlled PEVs can reduce the load change 
from one hour to the next in both cases. For the extreme values, 
the variation between different weather years is high (see Table 5). 

For the 2008 weather data, the highest ramp down is —14.41 GW 
versus —14.45 GW using controlled PEVs and the highest ramp up 
is 22.13 GW versus 17.90 GW. For the extreme ramp down value 
the used load management mechanism does not necessarily result 
in system benefits (see Table 5). 
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Fig. 12. Detail screen of the negative residual load when charging after the last trip and using demand-side management (DSM). 


D. Dallinger, M. Wietschel / Renewable and Sustainable Energy Reviews 16 (2012) 3370-3382 


-2500 


-5000 


-7500 


-10000 


Residual load (RS) change (RSq)-RS¢.1)) [MW] 


-12500 


-15000 


EZ Residual load (weather year 2006-2008) 


E Residual load + PEV using DSM (weather year 2006-2008) 


0 500 1000 1500 2000 2500 3000 3500 


3379 


—— 


4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 


Hours of the year 


Fig. 13. Sorted hourly load change of the residual load versus the residual load change using demand-side management (DSM) with plug-in electric vehicles. 


7. Conclusions 


In the presented paper, we investigated the capability of plug-in 
electric vehicles (PEVs) to balance intermittent renewable energy 
sources (RES) in a 2030 case study for Germany. The conclusions in 
detail are: 

PEVs provide a very high power/energy ratio: Compared with 
other storage devices, PEVs are able to offer a high total connec- 
tion power. A fleet of PEVs provide power totalling 54.12 GW (2030 
scenario with 12 million PEVs) with a relatively low usable amount 
of battery storage of 123.95 GWh (ratio 0.44). By comparison, Ger- 
man pumped-storage plants provide 7.76 GW with a rated volume 
of 224.31 GWh (ratio 0.035). 

Driving behavior restricts the use of mobile storage: The main pur- 
pose of PEVs is to fulfill mobility needs at equivalent costs to those 
of conventional vehicles. A cost-sensitive consumer will maximize 
the distance driven electrically in order to recoup the higher ini- 
tial investment. This implies high utilization of the battery and 
therefore reduces the time period available for load shifting and/or 
vehicle-to-grid services. PEVs are therefore utilizable as a short- 
time storage option (1-2 days) with limitations (e.g. infrastructure 
or consumer needs) on the load management time during the day. 

Price-based load shifting mechanisms require a feedback loop: 
The load management mechanism presented uses dynamic prices 
as a control signal for PEVs and takes the effect of intermit- 
tent renewable energy sources, mobility behavior as well as local 
grid restrictions into account. To avoid avalanche effects of auto- 
mated control a feedback of transformer utilization is included. 
This implies that PEVs have access to information about the reac- 
tion of other PEVs in the same distribution network. In practice, 
this could be realized in the form of distributed grid monitor- 
ing carried out by PEVs and/or the presented dynamic grid fee. 
Such a system does not yet exist but seems to be technically 
feasible. 

The consumer reaction on price signals is unclear: The economic 
incentives from electricity markets are low. The current base peak 
spread at the European Energy Exchange (EEX) is only about 
3 ct./kWh. In the presented approach, we assumed an infinite price 
elasticity of the consumer. Consequently, the PEVs’ load is able to 
be shifted to different time periods even in the case of very low 


incentives. Real consumer reaction in the case of the used load 
management mechanism is unclear. 

PEVs can contribute to balancing intermittent RES: To analyze 
PEVs’ contribution to balancing intermittent RES, three evaluation 
parameters were defined: the change in the minimum residual 
load, the percentage of negative residual load that can be con- 
sumed using load shifting, and the residual load change. For all 
three aspects, PEVs make a positive contribution to improving the 
integration of intermittent RES into the electricity system. 

The characteristic of RES intermittency greatly influences the 
results: The results are influenced by the generation share and the 
profile of intermittent RES, e.g. wind or photovoltaic. Especially 
for the percentage of negative residual load that can be consumed 
we found high variation depending on the RES characteristic. The 
expected higher utilization of power from German offshore wind 
parks during night-time hours supports the positive contribution 
of PEVs to integrating intermittent RES. The amount of power 
generated by wind turbines obviously depends on local climatic 
conditions and varies from country to country. Therefore, the pre- 
sented results are only valid for a German 2030 case study and the 
calculated values cannot be universally applied. 
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Appendix A. Appendix 


See Tables 6-8. 


Table 6 

Probability of travel. 
Mon Weekdays Fri Sat Sun 
62.73% 65.86% 64.94% 54.99% 40.66% 
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Table 7 
Probability for starting a trip (sliding average). 

Time Mon Weekdays Fri Sat Sun 

0 0.04% 0.06% 0.09% 0.17% 0.08% 
1 0.03% 0.03% 0.07% 0.14% 0.05% 
2 0.03% 0.04% 0.13% 0.17% 0.07% 
3 0.02% 0.03% 0.12% 0.12% 0.04% 
4 0.02% 0.03% 0.12% 0.12% 0.04% 
5 0.00% 0.03% 0.09% 0.07% 0.03% 
6 0.00% 0.03% 0.09% 0.06% 0.03% 
7 0.00% 0.03% 0.03% 0.04% 0.02% 
8 0.03% 0.03% 0.04% 0.08% 0.02% 
9 0.03% 0.02% 0.03% 0.08% 0.01% 
10 0.03% 0.02% 0.05% 0.07% 0.05% 
11 0.03% 0.02% 0.03% 0.06% 0.05% 
12 0.03% 0.03% 0.05% 0.07% 0.11% 
13 0.01% 0.03% 0.03% 0.03% 0.10% 
14 0.02% 0.03% 0.09% 0.05% 0.12% 
15 0.02% 0.03% 0.12% 0.05% 0.08% 
16 0.02% 0.04% 0.13% 0.05% 0.08% 
17 0.04% 0.06% 0.12% 0.03% 0.03% 
18 0.10% 0.15% 0.17% 0.02% 0.09% 
19 0.14% 0.18% 0.18% 0.04% 0.10% 
20 0.28% 0.30% 0.28% 0.08% 0.14% 
21 0.33% 0.36% 0.35% 0.11% 0.14% 
22 0.56% 0.53% 0.50% 0.19% 0.15% 
23 0.64% 0.59% 0.58% 0.21% 0.14% 
24 0.91% 0.88% 0.76% 0.25% 0.13% 
25 1.11% 1.03% 0.95% 0.30% 0.13% 
26 1.39% 1.40% 1.21% 0.39% 0.18% 
27 1.53% 1.54% 1.29% 0.37% 0.19% 
28 1.79% 1.97% 1.63% 0.49% 0.28% 
29 1.77% 1.93% 1.67% 0.57% 0.37% 
30 1.80% 2.02% 1.70% 0.78% 0.54% 
31 1.65% 1.77% 1.57% 0.90% 0.63% 
32 1.55% 1.78% 1.71% 1.31% 0.83% 
33 1.42% 1.46% 1.50% 1.42% 0.94% 
34 1.60% 1.57% 1.59% 2.07% 1.25% 
35 1.46% 1.41% 1.39% 2.10% 1.32% 
36 1.70% 1.52% 1.59% 2.69% 1.60% 
37 1.59% 1.38% 1.41% 2.54% 1.64% 
38 1.76% 1.53% 1.57% 3.27% 1.88% 
39 1.52% 1.30% 1.48% 2.93% 1.68% 
40 1.64% 1.41% 1.66% 3.24% 1.75% 
41 1.37% 1.30% 1.47% 2.81% 1.52% 
42 1.65% 1.53% 1.63% 3.25% 1.90% 
43 1.38% 1.31% 1.43% 2.70% 1.69% 
44 1.49% 1.45% 1.50% 2.76% 2.09% 
45 1.29% 1.37% 1.31% 2.43% 2.08% 
46 1.57% 1.68% 1.49% 2.76% 2.61% 
47 1.33% 1.53% 1.41% 2.34% 2.19% 
48 1.57% 1.71% 1.64% 2.51% 2.15% 
49 1.43% 1.60% 1.45% 2.20% 1.78% 
50 1.65% 1.79% 1.79% 2.46% 2.20% 
51 1.43% 1.49% 1.68% 1.93% 1.68% 
52 1.58% 1.53% 1.71% 2.07% 1.90% 
53 1.37% 1.34% 1.52% 1.63% 2.08% 
54 1.56% 1.65% 1.99% 1.91% 2.57% 
55 1.35% 1.40% 1.64% 1.46% 2.12% 
56 1.46% 1.49% 1.87% 1.77% 2.50% 
57 1.29% 1.43% 1.79% 1.49% 2.24% 
58 1.65% 1.70% 2.19% 1.96% 2.54% 
59 1.54% 1.38% 1.85% 1.61% 1.98% 
60 1.77% 1.64% 2.16% 1.78% 2.04% 
61 1.81% 1.58% 1.92% 1.50% 1.68% 
62 2.33% 2.00% 2.32% 1.68% 2.07% 
63 2.06% 1.84% 1.95% 1.17% 1.56% 
64 2.30% 2.26% 2.15% 1.26% 1.78% 
65 2.11% 2.12% 1.99% 1.05% 1.68% 
66 2.38% 2.68% 2.30% 1.11% 2.19% 
67 2.05% 2.27% 1.94% 0.87% 1.86% 
68 2.23% 2.40% 2.00% 1.06% 2.07% 
69 2.08% 2.12% 1.76% 0.94% 1.81% 
70 2.43% 2.42% 1.99% 1.25% 2.34% 
71 2.01% 1.85% 1.64% 1.15% 1.77% 
72 1.92% 1.92% 1.61% 1.36% 2.07% 
73 1.61% 1.71% 1.40% 1.13% 1.78% 
74 1.69% 1.69% 1.47% 1.41% 2.04% 


75 1.35% 1.27% 1.12% 1.06% 1.32% 
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Table 7 (Continued) 


Time Mon Weekdays Fri Sat Sun 
76 1.29% 1.26% 1.13% 1.19% 1.31% 
77 1.17% 1.15% 1.03% 1.01% 1.00% 
78 1.18% 1.15% 1.14% 1.10% 1.19% 
79 0.83% 0.92% 0.89% 0.76% 0.81% 
80 0.75% 0.84% 0.83% 0.77% 0.92% 
81 0.64% 0.69% 0.64% 0.52% 0.82% 
82 0.63% 0.66% 0.55% 0.46% 0.88% 
83 0.57% 0.50% 0.37% 0.32% 0.64% 
84 0.61% 0.52% 0.43% 0.34% 0.76% 
85 0.56% 0.46% 0.38% 0.34% 0.62% 
86 0.67% 0.57% 0.52% 0.40% 0.75% 
87 0.52% 0.44% 0.46% 0.32% 0.58% 
88 0.45% 0.44% 0.47% 0.41% 0.66% 
89 0.35% 0.34% 0.39% 0.45% 0.54% 
90 0.35% 0.35% 0.43% 0.51% 0.54% 
91 0.19% 0.20% 0.27% 0.42% 0.36% 
92 0.18% 0.18% 0.29% 0.43% 0.30% 
93 0.13% 0.11% 0.22% 0.33% 0.12% 
94 0.12% 0.11% 0.18% 0.24% 0.11% 
95 0.04% 0.06% 0.08% 0.12% 0.08% 
Table 8 
Probability of travel (cumulative). 
k Km Mon Weekdays Fri Sat Sun 
0 <2 0.22917949 0.21110176 0.224643 0.21110176 0.20020176 
1 <4 0.40713438 0.38936243 0.39172674 0.38936243 0.35961713 
2 <6 0.50709795 0.51071555 0.50926986 0,51071555 0.47829158 
3 <8 0.59882024 0.59218287 0.60156551 0.59218287 0.55480705 
4 <10 0.65561197 0.65367563 0.66401779 0.65367563 0.61194712 
5 <12.5 0.69032293 0.70196422 0.71391008 0.70196422 0.66313738 
6 <15 0.75472976 0.7635326 0.76389908 0.7635326 0.7251566 
7 <18.5 0.80461588 0.80302923 0.81311949 0.80302923 0.771473 
8 <20 0.83900508 0.83542293 0.84329914 0.83542293 0.81502836 
9 <25 0.88860334 0.88616926 0.90470785 0.88616926 0.86084511 
10 <30 0.92303152 0.92040934 0.92839391 0.92040934 0.90145597 
11 <35 0.94932388 0.94556082 0.94298257 0.94556082 0.91585863 
12 <40 0.96128432 0.9583068 0.96231069 0.9583068 0.92247889 
13 <45 0.96467203 0.96685348 0.96801052 0.96685348 0.92717667 
14 <50 0.9763209 0.97537706 0.972619 0.97537706 0.93510342 
15 <60 0.97966632 0.98305397 0.98302076 0.98305397 0.94832225 
16 <70 0.98367105 0.98741263 0.9866826 0.98741263 0.95747725 
17 <100 0.9908388 0.992556 0.99166396 0.992556 0.97390813 
18 <150 0.99442675 0.99554262 0.99609742 0.99554262 0.98236626 
19 <300 0.99876699 0.99794414 0.99809345 0.99794414 0.99491783 
20 >300 1 1 1 1 1 
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